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Abstract: We present formulas for the coefficients of 2-, 3-, 4- and 5-point master integrals for one- 
loop massive amplitudes. The coefficients are derived from unitarity cuts in D dimensions. The input 
parameters can be read off from any unitarity-cut integrand, as assembled from tree-level expressions, 
after simple algebraic manipulations. The formulas presented here are suitable for analytical as well as 
numerical evaluation. Their validity is confirmed in two known cases of helicity amplitudes contributing 
to gg — > gg, gg — > gH, where the masses of the Higgs and the fermion circulating in the loop are kept as 
free parameters. 
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B.3.31 The case k > 4 |4 



1. Introduction 

The unitarity method introduced in ||, ||| is designed to compute any scattering amplitude by matching 
its unitarity cuts onto the corresponding cuts of its expansion in a basis of master integrals Q with rational 
coefficients. Each of these coefficients can be determined quantitatively from prior knowledge of the master 
integrals and the singularity structure of the amplitude. 

As the master integrals form a basis for amplitudes, so the unitarity cuts of master integrals have 
uniquely identifiable analytic properties, and can be used as a basis for the cuts of any amplitude. There- 
fore, the coefficients of the linear combination can be extracted systematically through the phase-space 
integration (instead of complete loop integration). 

Recently, unitarity-based methods for one-loop amplitudes have been the subject of an intense inves- 
tigation, through different implementations of the cut-constraints @, ||. ^. [?]. || ||. [h| 11, [l2|, [l^, [D| [l6|, 



17],|8j, |19|, [2g,£l|,£2|. 

The holomorphic anomaly of unitarity cuts M, pi] simplifies the phase-space integration dramatically: 



cut-integrals can be done analytically by evaluating residues of a complex function in spinor variables [ 24 ] , 
reducing the problem of so-called tensor-reduction to one of algebraic manipulation. 

Accordingly, in [||, O], a systematic method was introduced to evaluate any finite four-dimensional 
unitarity cut, yielding compact expressions for the coefficients of the master integrals. This method was 
successfully applied to the final parts of the cut-constructible part of the six-gluon amplitude in QCD. The 
same method, based on the spinor-integration of the phase-space, was later extended for the evaluation of 



generalised cuts in D dimensions p^, 13, 14, 15], which is essential for the complete determination of any 



amplitude in dimensional regularization |2^, |26|, pTfl . 

In this paper, we carry out the extension to the massive case of the analytic results presented in 



[16], stemming from an original study of compact formulas for the coefficients of the master integrals [14]. 



Following the same logic as in [16], we now present general formulas for the coefficients of the master 
integrals which can be evaluated without performing any integration. These formulas depend on input 
variables (indices, momenta and associated spinors) that are specific to the initial cut-integrand, which is 
assembled from tree-level amplitudes. The value of a given coefficient is thus obtained simply by pattern- 
matching, that is by specializing the value of the input variables to be inserted in the general formulas. 
The implementation of the general formulas into automatic tools is straightforward, as done for the current 



investigation with the program S@M [28] 



In this paper, since the formulas for the coefficients are obtained via massive double cuts in D- 
dimension, we do not present results for the coefficients of cut-free functions like tadpoles and bubbles 
with massless external momentum (which can be expressed in terms of tadpoles as well). The coefficients 
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of such functions could be fixed either by imposing the expected UV-behaviour of the amplitude, as 
described in 26], or computed with other techniques applicable in massive calculations 0, U, H, H H • 
The paper is organized as follows. In section 2, we describe the structure of the decomposition of one- 
loop amplitude in terms of master integrals. In Section 3 we explain the double-cut integration with spinor 
variables, which leads to the formulas of the coefficients of the master integrals, presented in Section 4. In 
Sections 5 and 6, we apply our formulas to two examples of one-loop scattering amplitudes, respectively 
gH — * gg and gg — > gg, where the Higgs mass and the mass of the internal fermion (in both cases) are kept 
as free parameters. In Section 7, we present both analytical and numerical methods to obtain, finally, the 
explicit coefficients of the dimensionally shifted master integrals. In Appendix A, we record the translation 
between our basis of integrals and the ones used in the literature for the examples discussed in Sections 5 
and 6. In Appendix B, we present a proof of the decomposition into the dimensionally shifted basis, with 
rational coefficients independent of e. In other words, we prove that the coefficients given by our algebraic 
expressions will be polynomial in our extra-dimensional variable u. As a byproduct, we have produced 
equivalent and simpler algebraic functions for the evaluation of coefficients. 



2. Decomposition in terms of master integrals 

We define the n-point scalar function with non-uniform masses as follows: 1 

d 4 ~ 2e p 1 



I n (Mi,M 2 ,mi, . . . ,m n _ 2 ) = 



(2vr)4-2e {p 2 _ M 2 )((p _ R) 2 _ M | } _ Pj) 



2 — m 2 ) 



(2.1) 



Giele, Kunszt and Melnikov [22| have given the decomposition of any one-loop amplitude in D dimensions 
in terms of master integrals, represented here pictorially. 
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(2.2) 



Here, with reference to [22]: i) we have absorbed the residual D-dependence of the coefficients in the 
definition of the master integrals; ii) for ease of notation, we have given as understood the sums on the 



x For ease of presentation, we are omitting the prefactor i( — l) n+1 (47r) £,//2 (which was included for example in [B6| 
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partition of the n-points of the amplitude in the number of points corresponding to each master integral. 
Thus, the coefficients e, d, c, b, a in Eq. fl2.2|) are independent of D. 



If, on both sides of Eq.(2.2), we apply the standard decomposition of the D = 4 — 2e dimensional loop 
variable, L, in a four-dimensional component, t, and its (— 2e)-dimensional orthogonal complement, //, 



then the integration measure becomes 



d A ~ 2 *L = / <T'"ii I dH 



-2e. 



(2.3) 



(2.4) 



namely the composition of a four-dimensional integration and an integration over a (— 2e)-dimensional 
mass-like parameter. By taking the //-integral to be understood, the four-dimensional integration on both 
sides of Eq.(|2.2D, can be read as follows: 



V 




/ 



/ 



\ 




+ h(iJ 2 




+ a 



(o) 



i< 4 > \ 2 -5) 



where d n (p?) , c n (fj?) , and b n (/i 2 ) are polynomials of degree n in /i 2 , as discussed in Appendix | 

d 2 {n 2 ) = dW+dWf + d^{i?f , 

Cl ( M 2) =c (0) +c (2) At 2 ) 
&1 ( M 2) = 6 (0) + 5(2)^2 ^ 



(2.6) 
(2.7) 
(2.8) 



whereas tt(/j. 2 ) is non-polynomial in // 2 and corresponds to the coefficients of the reduction of the pentagon 
to boxes, which occurs in D = 4. 

The polynomial structure of d n (fj, 2 ),c n (fi 2 ), and b n (fi 2 ) is responsible for the dimensionally shifted 
integrals appearing in Eq.( |2.2D , because the //-integration can be performed trivially by absorbing the 



extra powers of fj, into the integration measure, according to [26]: 



/ 



'/ ''I' ,-.2.r 



(27T) 



-2c 



(^Y /(,/) = -6(1 - e )(2 - e ) ■ • • (r - 1 - 6)(4vr)' 



(27r) 2r -* 



(2.9) 



The presence of vr(// 2 ) in the coefficient of the four-dimensional box is a unique signature of the pentagon. 
We conclude that the reconstruction of the four-dimensional kernel of any one-loop amplitude, given in 
Eq.( |2.5| ), contains all the information for the complete reconstruction of the amplitude in D-dimensions, 
given in Eq.(f2~l]|), 

In the following pages, we present the general formulas of the coefficients of the box, triangle, 
7g 4 \ and bubble, \ obtained from the double cut of Eq.(|2 
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,(0) 



7r(// 2 ) + d 2 (M 2 



+ ci(/i : 




(2.10) 



Since the formulas for the coefficients are obtained via double cuts, we do not present the results for the 
coefficients of cut-free functions like tadpoles and bubbles with massless external momentum (which can be 
expressed in terms of tadpoles as well). Their coefficients could be fixed either by imposing the expected 



UV-behaviour of the amplitude, as described in [26], or computed with alternative techniques [17, [T^, U 



3. The double cut phase space integration 



In this section, we review the D-dimensional unitarity method [ P2j , 15] as applied in cases with arbitrary 
masses 14]. Our goal is to describe the structure of the cut integrand, from which we will directly read 
off the coefficients from the formulas in the following section. The formulas will be the massive analogs of 



the ones in [16|. 

Recall the phase space integration of a standard (double) cut in D = 4 — 2e dimensions. We use the 
usual decomposition of the D-dimensional loop variable, L, in a four-dimensional component, £, and a 
transverse (— 2e)-dimensional remnant, /x, 



L 



The integration measure becomes 

d 4 ~ 2 *L = 



d- 2e iJL I d 



T(-e) 



dfi 2 (fx 



2\-l-e 



d 1 



(3.1) 



(3.2) 



namely the composition of a four-dimensional integration and an integration over a (— 2e)-dimensional 
mass- like parameter. In order to write the four-dimensional part in terms of spinor variables associated to 
massless momentum, we proceed with the following change of variables: 



£ + zK, £ 2 = 0, 



(3.3) 



where £ is a massless momentum and K is the momentum across the cut, fixed by the kinematics. Accord- 
ingly, the four-dimensional integral measure becomes 



d^£ 



dz d 4 £ 5 + {£ 2 ){2£-K) . 



(3.4) 



The Lorentz-invariant phase-space (LIPS) of a double cut in the i^ 2 -channel is defined by the presence of 
two (5-functions imposing the cut conditions: 



J d 4 " 2 ^ = J d 4 ~ 2t L 5{L 2 - Ml) 5{{L -K)- M%\ 



(3.5) 



- 5 - 



Here Mi and M 2 are the masses of the cut lines. By using the decomposition of the loop variable in 
Eq. (|2,3| ), the four-dimensional integral can be separated, so that 

J d 4 "** = J dfi 2 (/ ,2 r l-e J d% m 

where the four-dimensional LIPS is 

d 4 <p = [ d 4 £5{£ 2 - Ml - 1?) 5((I- K) 2 - Ml - i?) . (3.7) 



The change of variables in Eq.(|3,3|), and the z- integration (trivialized by the presence of 5's), yield the 
four-dimensional LIPS to appear as 

J d A <j) = J d A £ 5 + (£ 2 ) i5((l - 2z)K 2 — 2£ ■ K + M\ — M|), (3.8) 

where 



_ (K 2 + M 2 - M 2 ) - ^A[K 1 ,M 1 ,M 2 ] - 4^ 2 M 2 

with 

A[K, Mi,M 2 ) = (K 2 ) 2 + {Ml) 2 + (Mf) 2 - 2K 2 Ml - 2K 2 M$ - 2M X 2 M 2 2 . (3.10) 



We remark that the value of z in Eq. (|3.9| ) is frozen to be the proper root (K > 0) of the quadratic argument 
of S(z(l - z)K 2 + z{Ml — Mf) — M 2 — ji 2 ), coming from 5(£ 2 — Ml — fj 2 ). For later convenience, one can 
redefine the /i 2 -integral measure as 



where the relation between u and fj, is given by 

_ 4if 2 /i 2 2 f uA[K,M 1 ,M 2 ] \ 
U = A[K, Ml ,M 2 r " = { 4^ ) ■ (3 - 11} 

We observe that the domain of u, i.e., u € [0, 1], follows from the kinematical constraints, as discussed in 
0- 

Finally, after the above rearrangement, the Z)-dimensional Lorentz-invariant phase-space of a double 
cut in the K 2 -channel can be written in a suitable form, 

J d 4 ~ 2e <S> = X {e,K,Mi,M 2 ) du u" x ~ e J d 4 (p, (3.12) 



where 



( KM M\ (4vr) e f A[K,M u M 2 ] ^ 
X (e,K,M 1 ,M 2 ) = — — - I — 2 ) , (3.13) 
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and where cftcj) was given in Eq,( |3.8[) . By using the definition of u given in Eq.( |3.1l| ), we can write 

a - - u 

z = , (3.14) 

where 



K 2 + M 2 -M 2 ^A[K,M 1 ,M 2 ] 
a= — 2 , (3= — 2 . (3.15) 

Notice that when M\ = M 2 = we have a = (3 = 1, thus reproducing the massless case. A useful relation 
between z and u is the following: 

. M? - M 2 „ . 

( 1 ~ 2z ) + K 2 = PVT^i- (3.16) 

This relation will be used in Appendix [b] to prove that that the coefficients given in this paper are poly- 
nomials in u, or equivalently fi 2 . As discussed in the previous section, this feature is essential for the 
straightforward reconstruction of dimensionally-shifted master integrals. 

The main feature of a double-cut LIPS parametrized as in Eqs ([3.12|J3.8l) is that the kernel of the 
integration is represented by the four-dimensional integral. In fact, the u-integration (or equivalently, the 
/^-integration) , is simply responsible of the rise of shifted- dimension master integrals. Thus, our interest 
in the extraction of the coefficients of the master integrals from a four-dimensional massive double cut, see 
Eq.(2.10), translates in focusing the discussion only on the J d 4 (j). 

The D-dimensional double cut of any one-loop amplitude is, in general form, 

j4-2e^. /i tree 4tree _ „ / ts K/r n,r \ I J„. ..— 1— e / j4 1 Atree .. /itrcc 



A % [ ee x A% cc = x (e, K, M 1 ,M 2 ) du iT 1 ' 6 J d 4 A % [ ee x A% ee (3.17) 

where A l £ ee and A^ ee are the two tree- level amplitudes on the left and right side of the cut. As discussed 
above, the kernel of the integration is represented by the four-dimensional part, 

d A 4> A^ x A% ee . (3.18) 

We proceed from the formula ( |3.8|) by introducing spinor- variables according to p3f| , 

d A £ 5(£ 2 ) = [ (I di) \l dl] [ t dt (3.19) 



and performing the integral over t trivially, with the second delta function. The general expression of the 
double cut integral will then be 

dV A l £ ee x A% ee 

/ d 4 £5 + (£ 2 )5({l - 2z)K 2 -21-K + M 2 - M 2 ) ^-=-^ } - (3.20) 

J U.M - K i) 2 - ™>i - » 2 ) 
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d 4 £S + (£ 2 )5((l - 2z)K 2 -2£-K + Mf- M| 



Uj (aj\£ + zK\bj 



X\ l {K 2 + M 2 1 -m 2 -2{£ + zK).K l 



{£ d£) [£ d£] / tdt5((l - 2z)K 2 + t (£\K\£) + Ml - M 



\[j{z{aj\K\bj]+t{£\Pj\£\) 



Ui(K? + Ml - ml -2zK-Ki + t (£\Ki\£]) 
2£- P, with P = la) [61. 



Here we have used £ 2 = M\ + jj?. Notice that \ a\£\b 

After using the remaining delta function to perform the integral over t, we have 

in+k 



{£ d£) [£ d£] (1 - 2z) + 



m 2 -m 2 \ {k 2 y+^ n]z?mM 



k 2 j (£\ K \£] n+2 nti 



where 



Pi = 


\ aj )[bj 




Rj = 




-2z) + 


Qj = 




-2z) + 



M\ - M 2 
Ml - Ml 



3 K ' 



K 2 + Ml -m 2 - 2zK ■ Kj 

K 2 J -3 ^ ^2 



Ki + 



K 



(3.21) 

(3.22) 
(3.23) 

(3.24) 



The vectors Pj,Rj,Qj do not depend on the loop variables, but rather only on the external kinematics, 



and especially on the momentum across the cut, K. Applying ( 3.21 ) to the master integrals, we find the 
following results. 



1) Bubble: k = 0,n + k = 0. 

(£ d£) [£ d£] [a- 2z ) + 

• (2) Triangle: k = 1, n + k = 0. 



Mf-M|\ (K 



K 2 



\K\ey 



(3.25) 



J (£ d£) [£ df] [ I - 2-.) + 
(3) Box: k = 2, n + k = 0. 



Ml - Ml 
K~ 2 



(£\K\£] (£\Q\£] 



(3.26) 



{£ d£) [£ d£] (1 - 2z) + 



Ml - M 2 
K 2 



(K 2 ) 



2\-l 



(£\Qi\£] (£\Q2\£] 



(3.27) 



These formulas are the extension to the massive case of the corresponding ones given in [|lq ]. We 
notice that the presence of the masses enters only the definitions of Pj , Rj , and Qj . Therefore the spinor 
integration performed in the massless case [|l(| is valid as well in this case. 



The expression of the cut-integrand in Eq.(3.21), with its indices, n and k, and its vectors Pj, Rj, 
and Qj is the key to constructing the coefficients. In the next section we present general formulas for the 
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coefficients of the master- integrals (boxes, triangles and bubbles), which depend on exactly these input 
parameters. Accordingly, given a specific amplitude (or integral), one can obtain its decomposition in 
terms of master-integrals without any integration. Every coefficient is obtained from the general formulas 
simply by substituting the input parameters characterizing the specific amplitude. These parameters are 



obtained by pattern- matching onto the reference form in Eq.( 3.21 ). 



4. Formulas for the coefficients of master integrals 



The coefficients of master integrals are obtained by the procedure described in the previous section, 
which is a straightforward generalization of the massless case ]l6|. We list the results in this section. 
In fact, the expressions take the same form as in the massless case; the mass dependence enters directly 
through the definitions (|3.23| ) and ( |3.24 ), and through these formulas into the definitions ( [4.2|) , (4.4). 



4.1 Box coefficient 



K s 




K - K r K a 

Figure 1: Double-cut of Box functions 



The formula for the coefficient of either of the box-functions with external kinematics as shown in Fig.^J is 



C[Q r ,Q s ,K] 



(K 



2\2+n 



^ \ (Psr,l\K\Psr,2] n+ Ylf-i^j {P s r,l \Qt\Psr,2 



+ {P S r,l <— ^ Psr,2 } , (4-1) 



where 



2/~>2 



A sr = (2Q S ■ Q r ) z - AQ Z S Q 



P, 



sr,2 



Q S + 



2Q 2 r 

-2Q S ■ Q r — \fAl 



2Q 



Qr 



(4.2) 



4.2 Triangle coefficient 



The formula for the coefficient of the triangle function with external kinematics as shown in Fig.y is 
(K 2 ) 1+n 1 1 



C[Q S ,K] 



2 (V^) n+1 (n + l)l(P s ,iPs,2 



n+l 
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K 




where 



Figure 2: Double-cut of a Triangle 



d n+l ( YfS(PB,l-TPs?\RjQ8\Pe,l-TP a t) 



dT n+i y U k =i t ^ {PsA _ T P Sj2 \Q t Q s \P sA - tP s , 2 ) 



(4.3) 



T = 



A s = (2Q S • Kf - AQ 2 S K 2 

PM = Q.+( - 2q ' 2 * +V5; i* 

P„ 2 = Qj + ( -^^-^ I A' 
Note that the triangle coefficient is present only when n > — 1. 



(4.4) 



4.3 Bubble coefficient 



K 




Figure 3: Double-cut of a Bubble 
The formula for the coefficient of the bubble function with the external momentum K, shown in Fig.||, is 

k n 



cm = (* 2 ) 1+n E^r^ <U(^EE(<- il) ( s )-t-f , ( s : 

q=0 \ r=l a=q 



=0 



where 



(4.5) 



(277 -If) 



|X-r?7l?7] 



, (4.6) 



r=0 



<f ;1) (^)- 



("I) 



b+1 



6! 



(Pr,l Pr,2 



d?_ ( 1 (P r ,i-rP r , 2 |77|P r ,ij 
>dr b \(t + l) (P r}1 - T P r ,2\K\P r>1 ] 



t+i 
~t+T 



(P r ,i - rP r , 2 |Q r 7;|P r ,i - TP r:2 ) b U]=i (Pr,i ~ rP r , 2 \Rj(K + s V )\P r 

{P Tjl - TP r>2 \ V K\P rtl - TP T , 2 ) n+1 U k p=l, P ^r (P,l - rP r ,2\Qp{K 



- rP r 



r,2) 

+ sri)\P r>1 - tP T)2 ) 



T=0 



(4.7) 
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H L i 1 




lm 
c 4 



+ c [12|3|# 




+ c [i\2\-m] 




+ C [12|3H] 



3 La 2 



Figure 4: Double-cut in the si2-channel for A(l + , 2 + , 3 + , H). 



\b+l 



d b ( 1 (Pr,2 ~ TP r ,l\v\Pr,2] 



t+1 



t+1 



(P r ,2 ~ rP r<1 \Q r r,\P r<2 ~ rP r ,i) b U%i (fi,2 - rPr^RjjK + sr^P^ - rP r<1 ) 

(Pr,2 - rP rtl \ V K\P r:2 - TP rA ) n+1 Up=l,p^r ( P r,2 ~ TP r ,l\Qp(K + S V )\P r>2 ~ tP t ,_ 



(4- 



T=0 



where A r , P r) i, P r ,2 are given by ( |4.4[ ), and r/, r/ are arbitrary, generically chosen null vectors. Note that the 
bubble coefficient exists only when n > 0. 



5. Example I: s 12 -channel cut of A(l+, 2+, 3+, H) 

In this section as well as the next, we check our formulas by reconstructing some helicity amplitudes 
contributing to gH — ► gg and — ► (75 at NLO in QCD, both known in the literature p9| , |26[] , We present 
our calculations in detail. 

Our first example is the si2-channel cut of A(l + , 2 + , 3 + , H). This amplitude was first computed in 
1 29]. Here, to facilitate comparison, we follow the setup of |27|, where the amplitude was rederived using 
unitarity cuts. At one loop, every Feynman diagram has a massive quark circulating in the loop. The 
quark mass is denoted by m. 

The s-channel cut of A(l + , 2 + , 3 + , H) admits a decomposition in terms of cuts of master integrals as 



shown in Fig. |4[ Its expression, given in Eq.(4.20) of [27|, reads 



^(l - * - ; 2 , 3~*~, H) \ s _ channel 



1 m 



St 



^%4%(m 2 + M 2 ) - s] 



(4vr) 2 - e v 2 (12)(23)(31) [2t(s - m% 

-^iT^ J 3 2) [4(- 2 + -s]- i/ 4 [4(m 2 + M 2 ) " s] .} (5.1) 



where s = s% 2 , t = s 2 3 , u = m 2 H — s — t. One can thus read the following values for the coefficients: 



c\ m = -co i(4(m 2 + ^ 2 )- S ) 



C[12|3|fl] 

C [1\2\3H] = c 
c [12|3ff] = 



1st 

M 

2t(s - m^) 



4(m 2 + ^ 2 ) - s) , 



4(m 2 + n 2 ) - s) , 



(5.2) 
(5.3) 

(5.4) 
(5.5) 
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with 



m 2 st 



°° (4vr) 2 - e v 2 (12)(23)(31) " (5 ' 6) 

5.1 The reconstruction of the coefficients 

We now show how to reconstruct the coefficients given above with our formulas from Section We fol- 
low the definition of the integrand given by |p7fl . By sewing the tree-level amplitude A^ ee (—L\, 1 + , 2 + , — L 2 ) 



and A i £ ee (—Li,3 + , H, —L2) given in Eq.(4.1-4.2) of [27], and using the Dirac equation for massive fermion 



it is shown that the four-dimensional integrand of the s-cut, C12, can be written as: 2 

C12 = cq,i— + 00,2777^ (5.7) 

where 

mK 2 S23 m[l 2] 

C °> 1 = ~v(K2-m 2 H )(12) (2 3) (3 1) C °' 2 = y/2v (1 2) (5 ' 8) 

Nx = m [4(m 2 + fi 2 ) - s 12 ] (5.9) 

N 2 = 8m(m 2 + /x 2 )^ • e+ + k A ■ e+) + v^m ^ffl (5.10) 

(1 3) 

D 2 = (4- kxf-^-m 2 , D 4 = (h + Uf - fi 2 - m 2 (5.11) 

We need to classify the contribution of N\ and N2 to each coefficient. Observe that N\ is independent of 
the loop momentum variable, so we consider it as a single term, 

N\ = m [4(m 2 + ^ 2 ) - s 12 ] , (5.12) 

while N2 is treated as three separate terms, 

N 2 = N 2 ,i + N 2 ,2 + iV 2>3 (5.13) 



where 



(1 1 2| 31 

JV 2 ,l = 8m(m 2 + /x 2 )(fc 4 • e£) = 8m(m 2 + // 2 ) V ' 1 J (5.14) 

v2(13) 



N 2 , 2 = 8m(m 2 + M 2 )(^i ■ 4) = -8m(m 2 + fi 2 )^^- (5.15) 

v2(13) 

/- (1|4|4|4|3] 

iV 2 ,3 = V2m — r . (5.16) 



Note that here we use "twistor" sign convention for the antiholomorphic spinor product, which is the opposite of the 
"QCD" convention followed by [El] [a; ?,] Rozowsky = -[xy] BFM . 
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By pattern-matching onto the reference form in Eq.( 3.2l| ), each integrand can be characterised by the 
parameters given in the following table. 



integrand 


n 


k 


?i = \Pi.m\ 


N l /D 2 


-1 


1 




N 2 ,i/(D 2 D A ) 


-2 


2 




N 2 ,2/(D 2 D i ) 


-1 


2 


|1> [3| 


N 2 , 3 /(D 2 D A ) 


-1 


2 


k 4 \3] {l\k 4 



(5.17) 



These data are the input values that we need in evaluating the formulas of the coefficients of the master 
integrals. 

From this table we draw the following conclusions. Since N\/D 2 has k = 1 and n = — 1, it contributes 
only to a triangle coefficient; whereas N 2 ^ {i = 1, 2, 3), having n < —1, contributes to both box and triangle 
coefficients. There are no bubble contributions at all. Thus we have already reproduced the absence of 
bubbles, Eq.( |5.5| ), without any calculation. 

To apply our formulas, we need to identify the definitions of K, K\, and K 2 . By inspection of D 2 and 
-D4, along with the fact of working in the s-channel cut, we choose the following consistent definitions: 



K = ki + k 2 ; Ki = ki; K 2 = -/c 4 . 
Since there is a single massive quark circulating in the loop, we have 

M\ = M 2 = rrij = m, 

and 



1 - y/T^uJl 



Am? 
S12 



From ( |5.18| ), we use the definition ( [3.24D to construct 

Qi = -(1 - z)ki - zk 2 

(l-z)^f K 

Using ([4.4D , we also set up the following quantities useful for triangle coefficients: 
Ai = {l-2zf{K 2 ) 2 
Pi,! = (1 - 2z)k 2 , Pi, 2 = -(1 - 2z)h 
A 2 = (l-2z) 2 (K 2 -m 2 H ) 2 

P 2A = -(l-2z)k 3 , P 2i2 = 



l-2z)^k 3 + (l-2z)(l-^)k 4 



in i 



(5.18) 
(5.19) 

(5.20) 

(5.21) 
(5.22) 

(5.23) 
(5.24) 
(5.25) 

(5.26) 
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5.2 The box coefficient c\ m 



The box coefficient c\ m takes contributions from N21, N 22 , ^2,3'- 
c\ m = cq,2 8m(m 



+co,2^C[Q l ,Q 2 ,K]^) (5.27) 



where C[Q r ,Q s , K], defined in Eq.(fO|), is 



(K 2 ) 2+ " ( Viitl^srAHPsrd 

2 ^(P sri ll^l^2]" +2 nt = l )t ^ s {PsrAQt\Psr,l\ 



C[Q r , Q S ,K} = ( - — -ii^-g ■" _ J — — ■ + {P sr>1 „ P sr<2 } ) . (5.28) 



. C[Q 1 ,Q 2 ,K]^ 

This term, corresponding to n = —2 is trivial, since k = 2 and A^,i has no dependence on the loop 
variable, 

C[Q l ,Q 2 ,K]W = i (5.29) 

• C[Q!,Q 2 ,K]W and C[Q U Q 2 , K}^ 

C[Q 1 ,Q 2 ,K] ( - 2 ^ and C[Q U Q 2 , K}^ both correspond to n = -l,k = 2. They differ only in the 
definition of P\. Therefore we can compute them in parallel, and specialize later to the corresponding P\. 
With n = — 1, k = 2, the expression is 

C[ Ql ,Q 2 , K] = ^l + {P*i,i " W) (5-30) 

^ V \-f~21,l| A l"21,2j / 

= " (1 " 2Z) ~ \{P^\K\P 21>2 ] + (P 21 , 2 \K\P 21 J + Z{ ~ 2K • Pl) 

if 2 ( {P^fr\P^^ ,_ 0K . 

~ [ J 2 i (P21,l|^l^21,2](P21,2|^|P21,l] ) + { 1] 

For \P\) = |1) , \P\\ = |3], (so that, for any S, one has 2P\ ■ S = — (1\S\3\), one obtains 



C[Qi,Q 2 ,K] (2 ' 2) = 

ET 2 / (l-2z^s 23 (s 23 + K 2 -mj I )(l\2\3] \ ( {\-2zf 2 2 A" 1 

= - (1 - 2z) ^r^o^r JUi^r 23(S23+ ^ - mH) ) +z{im 

= (5.31) 
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For |Pi) = # 4 |3] , \Pi) = j6 4 |l), ( so that, 2Pi • 5 = - (l|/c 4 Sfc 4 |3]), one gets 

c[Qi,g 2 ,K] (2 ' 3) = 

-m^ (5.32) 



<1|2|3] _ 2 



The result for c\ m 



The total coefficient of our box is: 

lm o / 2 2m + 8m(m 2 +/x 2 ) ,(1|2|3], y/2m . 2 <1|2|3] . 

c\ m = c ,2 8m(m 2 + // 2 )fc 4 • e+ - C0j2 — fii(U-U) + ^ i2 __(_ m ^lJ_U) 

m 2 if 2 s 2 3 



2u (1 2) (2 3) (3 1) 



[4(m 2 + - m^] (5.33) 



Multiplying by — i/(47r) 2 e , to account for the difference in the definitions of master integrals, we confirm 
the result of j27j]. 

5.3 The triangle coefficient C[i|2|3#] 

The coefficient cn|2|3fn gets contributions from N\, N 2i2 and ^2,3: 

C[1 |2|3H] = 0),! N, ClQuK^-co /™^ 2 ^ C[ Ql ,K]^ + co, 2 ^ C[Q lt K\™ (5.34) 

where the general triangle coefficient, given in Eq. fl4.3j ), reads 
(K 2 ) 1+n 1 1 



C[Q S ,K] 



A s )«+i (n + l)!(P S)1 P s , 2 ) n+1 

. (5.35) 



d n+1 ( Uj=i (PsA ~ tP s , 2 \R 3 Q s \P s + - rP s , 2 ) 

X TT T 2 V 1 Ps 1 *-* Ps 2 [ 



r=0 



<?[Qi,#] (1) 

We have already observed that the iVi term is trivial. Here is how that shows up in our formulas. 
Read C[Q S , K] for s = 1 and k = 1, n = —1., and no Rj. 
The term inside the parentheses degenerates to 1. 

CiQuKl® =i(l + l)| T=0 = l. (5.36) 
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. C[Q X ,K]W andC[Qi,A]( 2 > 3 ) 

As we said already, and A2,3 differ in the definition of P\. Therefore we start by manipulating 
the general formula, and only at the very end we specialize each contribution using the corresponding P\. 

Since n = —l,k = 2, there is no derivative at all, so we can set r = from the beginning: 

c\q k] = - ( iMMAM 1 iMMAM ) (5 37) 

2\U k t=1 ^ s (Ps,l\QtQs\Ps,l) Ut=l^s(Ps,2\QtQs\Ps,2))' 

The one-mass triangle (l|2|3i?) corresponds to the value s = 1. 



C[Qi,K] - -- [ 1 z m + 1 j m ) (5-38) 



For |Pi) = |1) , \Pi] = |3], one gets 



CiQuKfV = lUffl (5.39) 
2 s 2 3 



For \Pi) = ^4 1 3] , \Pi] = ^4(1), one obtains 

C\G Km- 1 A2|4[3](1|411] . (1|4|3](1[4|2] \_ / m 2 H \ 

c[QuK] A (2| 4 |i] + vm )- {im v~^) (5 - 40) 

• The result for Ch^isa] 

The total coefficient of triangle (l|2|3ff) is: 

mK 2 s 2 3 x Al 2 2n 1 m[12] 8m(m 2 + fi 2 ) (1|2|3], 

C W 2 ^ = 2) (2 3) (3 I)" t 4(m + ^ } " ^ " V2v(12) 72(13) 

m [1 2] \/2m / m 2 ^ 



<1|2|3] 1 



V2v(l 2) (1 3) N 1 1 1 V 2s 2 3 

m 2 K 2 (2s 23 + AT 2 - m^) 
"2w(A: 2 -m^)(l 2) (2 3) (3 1) 



[4(m 2 + n 2 ) - m 2 H ] (5.41) 



Multiplying by i/(47r) 2 e , to account for the difference in the definitions of master integrals, we again 
confirm the result of |27 



5.4 The coefficient C[i2|3|#i 

The coefficient cr 12 |3|_ffi gets contributions from and A2,3, therefore can be written as 

cww = -co /^+p C[Q 2 ,K]^ +C0>2 f^ C[Q 2 ,K]™ ( 5 .42) 
• C[Q 2 ,K]^ and C[Q 2 ,K]W 
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1 L i 2 




Om 
c 4 



+ C 



[23]4|1] 




+ C 



[2|3|41] 




+ c [12|34] 




4 L 2 3 



Figure 5: Double-cut in the S23-channel for A(l ,2 ,3 + ,4 + ). 



The two-mass triangle (12|3|iT) corresponds to the value s = 2, and its coefficient can be obtained 
from Eq. 



C[Q 2 ,K\ = - 



2 V (3|1|2|3) [3|1|2|3] 
By using |Pi) = |1) , |Pj] = |3], one gets 

By using |Pi) = 13] , |Pi] = ^6 4 1 1) , one obtains 

LV ' J 2 V (3|1|2|3) [3|1|2|3] / V K 2 H 2s 23 



(5.43) 



(5.44) 



(5.45) 



The result of cn2|3|iT] 
The total coefficient of triangle (12|3|iJ) is: 



C[12\3\H] 



8m(m 2 + n 2 ) ( m\\ < 1 1 2 1 3] V2m ( m 2 H \ 2 (1 

m 2 (K 2 — m\ 



2|3] 



2S 23 



U4(m 2 + ,i 2 )-m|] 



(5.46) 



2t> (1 2) (2 3) (3 1) 

Multiplying by i/(4-7r) 2_e , to account for the difference in the definitions of master integrals, we again 
confirm the result of |]27| . 



6. Example II: ^-channel cut of A(l ,2 ,3+,4+) 

Our second example features a non- vanishing bubble coefficient. We study the t-channel cut of the 
gluon amplitude A{1~ , 2~ , 3 + , 4 + ), with a massive quark circulating in the loop. (As usual, t = S23.) 

The i-channel cut of A(l~ , 2 _ , 3 + , 4 + ) admits a decomposition in terms of cuts of master integrals as 
shown in Fig. [|, and its expression was given in Eq.(5.33) of |{26[j. After converting that expression into 
our basis of D-dimensional master integrals, as done in Appendix [A|, it reads 



A 



fermion 



(l-,2-,3+4 H 



(1 2) 2 [3 4] 2 (2 



t— cut 



St 



J 2 [l] + -Um l + //I - ~h[m 2 + A 
3 6t s 



2t 



+ -h[(m 2 + fj,y] - th[m 2 + fx 2 ] 
s 



(6.1) 



t— cut 
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One reads the following values for the coefficients: 



Om 

c 4 



C [23|4|l] 
c [2|3]41] 

c [23|41] 



r„ ( ^(m 2 + //) -t 1 (///" + fl 2 








co ( \ + ^-W + lx 2 ) --{m 2 + ll 2 ' 
3 3t s 



with 



co 



(1 2) 2 [3 4] ; 

7t 



(6.2) 

(6.3) 
(6.4) 

(6.5) 
(6.6) 



6.1 The reconstruction of the coefficients 



We now apply our formulas of Section |] to construct the coefficients given above. We follow the 
definition of the integrand given by ]26]| . By sewing the tree-level amplitude A^ ee (— L%, 2~, 3 + , L 2 ) and 
^4 4 ree (— L 2 , 4 + , 1~,L\) given in Eq.(2.3) of p6fl , and using the Dirac equation for a massive fermion, it is 
shown that the four-dimensional integrand of the i-cut, C23, can be written as: 3 

2iVi + N 2 



23 



D 1 D 2 



with 



m = - r (i\e 1 \4] 2 (2\e 1 \3y 



^23 



iV 2 
£>2 



«23 

{h + h 



(1 2) [3 4](1K 1 |4](2|£ 1 |3] 



(*i - ^ 2 ) 2 - ^ 2 



m 



(6.7) 

(6.8) 

(6.9) 

(6.10) 
(6.11) 



By pattern-matching onto the reference form in Eq.( 3.2l| ), each integrand is characterised by the 
parameters given in the following table. 



integrand 


n 


k 


\Pl)[Pl\ 


1^2) [P 2 | 


\Ps)[Ps\ 


1^4) 


N 1 /(D 1 D 2 ) 


2 


2 


|1)[4| 


|2)[3| 


|1)[4| 


|2)[3| 


N 2 /(D 1 D 2 ) 





2 


|1)[4| 


|2)[3| 







(6.12) 



We define 



K = k 2 + k 3 ; Ki = - 
P\ = P3 = A1A4, P 2 



= p 4 



k 2 -- 
A2A3. 



(6.13) 
(6.14) 



Recall that here we use "twistor" sign convention for the antiholomorphic spinor product, which is the opposite of the 



"QCD" convention followed by 



[^Bern-Morgan = _ r^, y] BFM , 
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Moreover, since we have a quark of mass m circulating in the loop, we take 

Mi = M 2 = mj = m. (6.15) 

Then, by applying (|3.9| ), we find 

m 2 + fi 2 

z(l-z)= R2 . (6.16) 

For the N\ term, n = 2. For the N 2 term, n = 0. Both terms give boxes, triangles and bubbles. 
From the definitions ( |jD , (|3.24[) we have 

Qi = (1 - z)fci + zfc 4 , Q 2 = -(1 - z)k 2 - zk 3 (6.17) 
R 1= R 3 = -(l - 2z)AiA 4 , R 2 = R 4 = -(1 - 2z)A 2 A 3 . (6.18) 



Further, the quantities defined in (4.2), ( |4.4j ) become 

A 12 = (1 - 2zfs l2 {K 2 + s 12 ) - (1 - 2z) 2 K 2 s 12 (6.19) 

Ai = (1 - 2z) 2 (K 2 ) 2 , (6.20) 

P 1( i = -(1 - 2z)h, P 1)2 = (l-2z)fci, (6.21) 

A 2 = (l-2z) 2 (iT 2 ) 2 (6.22) 

P 2 ,i = (1 - 2z)fc 3 , P 2 , 2 = -(1 - 2z)A: 2 . (6.23) 

6.2 The box coefficient c^ m 

The box coefficient C4 receives contributions from both iVi and N 2 , and can be correspondingly 
decomposed as: 



,0m _ 2 



-C[Q 1 ,Q 2 ,^]W+ (1 2 ^ 34] C[Qi,Q 2 ,^] (2) • (6-24) 

We discuss the computation of C[Q\,Q 2 , K]^ and C[<5i, Q 2 , if]*- 2 ) in detail, starting from the expression 
given in Eq.([4.1|). 

• C[Q U Q 2 ,K}W 

For the Ni term, with n = 2, the expression is given by 



(K 2 ) 4 ( (P 2 i,i|Pi|P 2 i, 2 ] 2 (P 2 i,i|P 2 |P 2 i, 2 ] 2 (P 2 i, 2 |Pi|P 2 i,i] 2 (P 2 i, 2 |P 2 |P 2 i,i] 2 ^ 



For analytic simplification, the following trace identity is helpful. 

{Px\R\P 2 ] (P 2 \S\Pi] = Tr (^—^ /Pi ft /P 2 JS] = ^(PxPPaS) 
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In terms of vectors, 

tv^{V 1 V 2 V 3 V 4 ) = i((2Vi • V 2 )(2V 3 • V 4 ) + {2V l ■ V 4 )(2V 2 ■ V 3 ) - (2V l ■ V 3 )(2V 2 ■ V 4 ) - Ue^V^V^Vf) 
The coefficient can then be expressed in terms of traces, and evaluated as follows. 
C[Q 1 ,Q 2 ,K]W 

_ (K 2 ) 4 ((tr_(P 21 , 2J Fs:P 2l! ii? 1 )tr_ 2 (P 21)2J Fs:P 2l! iii! 2 )) 2 + (tv^P^KP^R^ tvJ(P 21A KP 21:2 R 2 )) 2 ) 



2(tr _(P 21A KP 21t2 K)Y 

(K 2 ) 4 z 2 (l-z) 2 [3 4] 2 (1 2) 2 



(6.26) 



'12 



• C[Q U Q 2 ,K)W 

For the N 2 term with n = the expression is given by 

C[ Ql ,Q 2 , *]<"> = (Jfl ( (^Ml^^Mlg^ + { p 21j1 „ p 212} > ) . (6 . 27) 

z \ ("21,1 1 K21,2j / 

Combining the two terms over a common denominator, we have 

C[Q U Q 2 ,K]W 

(K 2 ) 2 (tr_ (P 2 i, 2 ^P 2 i,i^i) tr_ (P 2 i, 2 KP 2hl R 2 ) + tr_ ^21,1*^21,2^1) tr_ (P 2 i,iKP 21 , 2 R 2 )) 

2(tT_(P 21A KP 21 , 2 K)) 2 

= (j^m^M^ (6 . 28) 

• The result of C4" 1 

We add our two contributions together and replace z using ( |6.16| ). The total box coefficient is thus 

0m _ 2 (if 2 )V(l - z) 2 [3 4] 2 (1 2) 2 (1 2) [3 4] (if 2 ) 2 z(l - z) [3 4] (1 2) 
° 4 " (K 2 ) 2 s 2 2 + K 2 s X2 (b ' 2yj 

= (m 2 +/, 2 )[3 4] 2 (f 2) 2 / _ 2{m 2 +^ 2 ) \ ^ 
su \ s 12 J 

6.3 The triangle coefficients cr^um and C[2|3|4i] 

Both terms exhibit the symmetry of the amplitude, so our two triangles are not independent. 
The triangle coefficients cr 2 3|4|i] and cr 2 |3|4i] receive contributions from both N\ and N 2 , and can be 
correspondingly decomposed as: 

c [2 3|4|i] = -j^y 2 C[Qi, K] W + M^. C [Q lt Kp (6.31) 
c [2|3|4i] = -j^yC[Q 2 ,K]V + {12) K f 4] C[Q 2 ,K]M (6.32) 
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We discuss in parallel, first the contribution due to Ni to both coefficients, namely C[Q\,K]^ and 
C[Q2, K]^ 1 ' , and later the one due to iV 2 , namely C[Q±, K]^ and C[Q2, K]( 2 ' , where the triangle coefficient 
was given in Eq.([4.3|), 

• C[Q 1} K]W and C[Q 2 ,K]W 

Since the N\ term with n = 2, the triangle coefficient expression is given by 



C[Q l5 iT]« 



(K 



2\3 



1 1 d 3 / Tl J 4 = i<4-rl|iZ J -Qi|4-rl) | UU (1 ~ t4\RjQi\^ ~ ^) 

(l-r4|Q 2 Qi|l-r4) 



2 (VAT) 3 3! (4 l) 3 dr 3 I (4 - rl|Q 2 Q x |4 - rl) 



+ 



(1-2*) 1 d 3 /[4|g 1 |4) 2 (4-rl,2) 2 [3|g 1 |4-rl) 2 , t 4 [4|Q 1 |4) 2 (1-t4,2) 2 [3|Q 1 |1-t4) 



2 3! (4 1) dr 3 



(4-rl|Q 2 0i|4-rl) 



+ 



(1-t4|Q 2 Qi|1-t4) 



T = 



T=0 



(1 - 2z)(l - z) 2 (if 2 ) 2 d 3 / (4 - rl, 2) 2 [3|Qi|4 - rl) s 



12(4 1) 



dr 3 I (4-rl|Q 2 Q 1 |4-rl) 



T=0 



(l-z) 2 (if 2 ) 2 d 3 / ((4 2) - r (1 2)) 2 ((1 - z) [3 1] + tz [3 4]) 2 
12 dr 3 V -(l-z)(4|3|l] +TS 3 4 + r 2 2(l|3|4] 



r=0 



(6.33) 



A similar calculation shows that 



= 0, 



(K 



2\3 



1 1 d 3 / n j 4 =i(3-r2|^Q 2 |3-r2) ]Jj=i (2 - r3|^Q 2 |2 - r3) 

/A^) 3 3! (3 2) 3 dr 3 1 (3 - r2|Q 1 Q 2 |3 - r2) + (2 - r3|QiQ 2 |2 - r3> 



/ r=0 

(6.34) 



which can also be seen by the symmetry of the amplitude and the cut. 

• C[Q!,K]W and C[Q 2 ,K}^ 

For the iV 2 term with n = 0, the expression is simpler as 

C[Q U K]M = (i-^Xl-^ 2 d f (4-rl 2) [3|Qi|4-rl) + r 2 (1 - r4 2) [3|Qi|l - r4> 



2 dr V (4-rl|Q 2 Qi|4-rl) (1 - t4\Q 2 Qx\1 - t4) 

(l-2z)(l-z)JY 2 d /((4 2)-r(12))([3|Q 1 |4)-r[3|Q 1 |l))' 



T = 



dr 



(4-rl|Q 2 Qi|4-rl) 



r=0 



(6.35) 



A similar calculation shows 



r\n K ^- {1 ~ 2z) d f (3-T2 1)[4|Q 2 |3-r2)[3|Q 2 |3) r 2 (2 - r3 1) [4\Q 2 \2 - r3) [3|Q 2 |3) \ 
J 2 dr V (3-r2|Q!Q 2 |3-r2) + (2 - r3|QiQ 2 |2 - r3) J 



0, 



T=0 

(6.36) 



which can also be seen by symmetry. 
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The results of c^um and cr 2 |3|4i] 
Every term vanishes separately, so 



'[23|4|1] = 0, C[2|3|41] = 0. (6.37) 



The vanishing results for triangle coefficients is not obvious from the beginning. We suspect that there 
should be a more directly physical argument to see this point. 

6.4 The bubble coefficient cr 2 3|4i] 

The bubble coefficient cr 2 3|4i] receives contributions from both N\ and iV^, and can be correspondingly 
decomposed as: 

-[23,41] = -^[K]^ + { -^^ClK]V ( 6 .38) 

There is one subtlety regarding the calculation of the bubble coefficient. The formulas involve an 
arbitrarily chosen, generic auxiliary null vector r/. If r\ coincides with one of the Ki, we need to use a 



modified formula, given in Appendix B.3.1 of [16]. In this example, we illustrate both options. First, we 
show the result with a generic choice of rj; second, we use the formulas for the case rj = K\. Both are 
suitable for numerical evaluation, while the special choice of rj may simplify the analytic expression. We 



will find that the two results agree with each other, as well as with 26] 



6.4.1 Generic reference momentum 7] 

Let us start with the formulas for generic rj, given in Eq,(|4.5D, There are two terms we need to calculate. 



• C[K]W 

For the first term, with N\ in the numerator, and n = 2, the coefficient is 

q=0 y ' \ r=l a=q 



, (6.39) 

s=0 



where 



d 2 ( i (2r? -Kf-i n*=i<w^+^)K> 



2 > 2 -* y ' dr* \2[v\mv}HZ-q)(K 2 )^(£ v ) 3 nl =1 (£\Q p (K + SV )\£) 



\\e)-*\K-rrj\ V ] 



, (6.40) 

r=0 



■;a- q ;l) (s) = (~l) a ' g+1 d°^_ ( 1 (Pr,l-rP r ,2MPr,lf- C 

2 ~ a (a - q)\JK r a - q+X (P r>1 P r , 2 > a ^ dT-i {(3 -a) (p r>1 - T P r , 2 \K\P r>l f- 

(P r ,i - TP r , 2 \Q r ri\Pr,i ~ rPr,2} a ' q U 4 j= i ( P r,l ~ rP r ,2\Rj(K + sr,)\P r ^ - rfi, 2 ) \ 
(P rA - TP ri2 \ V K\P rA - TP r , 2 f Ul=i, P ^r (Pr,l - rP r ,2\Q P (K + ^)|^r,l - rP r , 2 ) ) 



-22- 



d a ~ q ( 1 {Pr,2-TP r ,M P r,2] 



3— a 



(a - g)!^A; a ^ +1 (Pri Pr,2) a ~ Q dra ~ q V ( 3 " «) <^,2 - rP,,i |if|P r , 2 ] 



3— a 



<fi,2 ~ rP r ,i|Q r 77|P r , 2 - rP rA ) a - q U 4 j= i (fi,2 - rP r ,i|^(K + sr?)|P r , 2 - rfi.i) \ 

X (Pr,2 ~ TP r>1 \ V K\P r>2 - TP rA f Ul=l^r <^,2 " tP T)1 \Q p (K + S V )\P r>2 - T P r>1 ) J 

After making some substitutions, and considering the summation ranges of a and q, we get 



T=0 



(6.42) 



(6.43) 



-1) 



a—q+l 



d a-q 



3— a 



1 <l-r4|7/|l] 

3— a 



(a-(?)!VAT a ^ +1 (4 l) a -'?dr a -'? ^(3 -a) (l- r 4|K|l] 



(1 - TA\Qm\l - r4) a ~ g {-tA\Ri{K + - t4) z (1 - r4| J R 2 (K + srj)\l - t4)^ 
(1 - tA\t]K\1 - r4) 3 (1 - t4\Q 2 {K + s?7)|l - r4) 



T=0 



(6.44) 



g(2;a-g;l)/ •> _ 

2 ' 2 ~ a (a - g)! VA^ a_9+1 (3 2) a - q dr°-9 \ (3 - a) (3 - r2|^|3] 3 - a 

, (3 - t2|Q 2 ??|3 - t2) q ~' ? (3 - r2| J R 1 (if + sr/)|3 - t2) 2 (3|fl 2 (ir + S7?)|3 - t2) 2 < 
(3 - T2|r/if|3 - r2) 3 (3 - r2|Qi(if + srj)\3 - r2) 



(-1) 



a-q+1 



1 <3-r2H3] 



3— a 



r=0 



(6.45) 



sg a _7 ;2) ( s ) = o. 



(6.46) 



C[if]( 2 ) 

For the second term iV 2 with n = 0, the expression is much simpler: 



where 



r=l 



/ (2rj • K) {t\RiK\t) (£\R 2 K\£) \ 

i^ 2 MtfQiM N^Q 2 |ry] 

K 2 [r? 1] [7/ 2] 

1 (4|t?|4] (4|fliif|4) (A\R 2 K\A) 
VA7(4|K|4] (4|^|4) (A\Q 2 K\4) 



s=0 



\l)^\K\rt] 



(6.47) 



(6.48) 
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1 M](4 2) 
K 2 [77 1] (4 3) 



<f 2) (* = o) = o. 



<6° ;1) ( S = o) 



<f 2) (* = o) = o- 



1 <3|tt|3] (3|i?iijf|3) (3|^ 2 K|3) 
"^(31^13] (317/^13X3^1^13) 
1 [77 3] [4 2] 



K 2 [rj 2] [1 2] 



(6.49) 
(6.50) 

(6.51) 
(6.52) 



• Results: 

We have used the numerical routines of S@M p8| to show that while each single term B entering 
Eq.(6.39) is //-dependent, their combination is indeed independent of the choice of 77. The choice 77 = k% is 



found to be convenient. (Note that k% is not proportional to either of K±, K2.) Therefore we set 77 = ks, 
and we obtain the following analytic result. 



3[12|34] 



2 r (1 3) [1 3] (k 2 ) 3 / 1 



+ 



(K 2 ) 2 \ [12] 2 (4 3> 2 \2 
K 2 [l 3] (1 3) / 1 2 



4z(l - 2) 



*(f> 4 a , (1 _ I) + <^) 4 



3[1 2f (4 3) z 



6 [1 2f (4 3) 



(4 3) 2 [12] 2 \6 3 



1 



(- - -z(l - *)) (1 3)^ [1 3]' + (- - 3z(l - z))K 2 (1 3) [1 3] + 



(1 2) 2 [3 4]' 
^H 2 



-4(m 2 + n 2 )s + 6(m 2 + /i 2 )i - 2st) 



(1 2) [3 4] [43] 
K 2 [1 2] 

(6.53) 



We have used the relations K 2 = t, (1 3) [1 3] = — s — t. 

6.4.2 Special choice of 77 

Alternatively, we discuss the calculation of the bubble coefficient by using the special choice of 77 = 
K\ = —k\ from the beginning. With this choice, we need to use formulas for the B which are slightly 
different from the ones used in the previous section. They are given in Appendix B.3.1 of [16|. 

Our convention for the spinors is: 



\v) 



N = -|i]. 



(6.54) 



c[K ]n = ± t}l£. U°i_ q{s) + (B^::f\ s ) - B^:f\s))) 

q=0 y ' \ r=2 a=q / 



(6.55) 



s=0 



jn+1 



[rj\r]K\rj] ' 



n-l 



dT n+i y (1 _ 2z ) - sz (t + i)(n + 1)! (i v ) n+2 U k p=2 (*\Qp{K + *n)W 



(6.56) 
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Since k = 2, we can directly set r = 2: 



s£ a;1) (*) 



i 



d a / [1 3] 



(n)»+ 1 - t -«[i4] n+2 (4<) n+2 J 



K)=|3)-r|2> 



t+1 



(1 - 2z) - sz (1 - 2z) a + 1 (K 2 ) a +*+ 2 a! (3 2) a dr« I + 



(6.57) 



s£ a;2) (*) - 



(-ir 



d a / [1 2] 



(£ 1) 



n+l—t—a ^ ^j«+2 ^ £} n+2 I 



K)=|2)-r|3> 



t+1 



(1 - 2z) -sz(l- 2z) a + 1 (i^ 2 ) a+t+2 a! (3 2) a dr a \ (t + 1) 



Now we proceed to evaluate. 



(6.58) 



For the N2 term with n = 0, the evaluation is simple. In particular, there are no derivatives in s, so 
we can set s = directly. 



s=0 



Choosing 77 = 3, we have 

<2(* = o) = 



[1 2] (3 4) 
so the total contribution comes to 



[43] 
[1 2] 



(6.59) 



1 , « ;1) (* = 0) = ^M|A < 2) (, = 0) = (6.60) 



(6.61) 



C[K\W 



For the A' - ! term with n = 2, the calculation is a bit more involved. 

q=0 q ' V a=q ) 



s=0 



For the various terms, we have 

1 d 3 /(l-2z) 4 (l + s) 2 



(3 - g)3! (3 4) 2 [1 2] 2 dr 3 \ (1 - 2z) 



(1 - r) 2 (i^ 2 (l + s) - t{K 2 + a [3 1] (1 3))) : 



(6.62) 



(6.63) 



((1 + a)(l - 2z)K 2 + r((zs + 2z - s - l)K 2 - s(l - 2z) [3 1] (1 3)) + r 2 s(l - z) [3 1] (1 3)) 
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s£^ ;1) to 



(1 + s) 2 (l - 2z) 3 - a+ i (-IT' 1 (3 2) 2 [1 3] 



i3— a 



(1 - 2z) - sz 



[K 2 f-i(a-q)\ 



(6.64) 



d a-q n(i _ z ) [! 2 ] + TZ [l 3]) a -"((l 3) - r (1 2)) 3 -<?(K 2 + s [3 1] (1 3) - sr [3 1] (1 2)) 



(3 -a) [1 4]" ((4 3>-r<4 2)Y 



The term B^_^' 2 \s) vanishes after taking the derivatives with respect to s. The reason is the following. 
Notice that 



B 



(2;a-<};2) 
n,t 



to 



(l + s) 2 (l-2z) 3 - a+< ? (-I) "' (3 2) 2 [1 2]* +1 d*-« ( t 2 (e\Q 2 \l] a - q [3\K + sr]\e) 



(1 - 2z) - SZ (K2)a-g+t+2( a _ q )\ (3 2 ) a - ( ? tfr°-9 ^ (t + 1) (£ l}- 1 "*-^? [1 4 ] 2 (4 

We can see that the r-derivative vanishes unless a — q = 2, in which case we get 

(l + s) 2 (l-2z) (3 2) 2 [1 2] m / 1 (-z) 2 (2|3|l]V[3|l|2) 2 \ 



(6.65) 



K>=|2)-r|3> 



(1 - 2z) - sz (K2)4+i (3 2 ) 2 I it + 1) (2 I)- 3 "* [1 4] 2 (4 2)' 



(6.66) 



However, the condition a — q = 2 implies a = 2,q = 0. Therefore we can set s = 0, and the expression 
vanishes: 



eg 2;2) to = o. 

Now we collect the results of ( |6.63| ) , ( |6.64 ) and ( |6.67D . We take rj = k% . Define 

[1 3] (1 3) 



Ci 



[K 2 ) 2 [l 2] 2 (4 3) 2 ' 



Let us begin with the terms with q = 0: 
K 2 (\-2z) 2 



B®(s = 0) 



B. 



(2;0;1) 
2.2 



3(3 4) 2 [2 l] 2 
1 



(6.67) 



(6.68) 



(6.69) 



( s = ) = d(--)(l-2z) 2 [1 3] 2 (1 3)' 



00 1;1) ( S -O) (', ( --(1 -2:r[l .">|-l :!:-!-(----. + :tf)h"[\ '■'>] [ I :$, 

? (2;2;1) 



1 



,3 9 



2 V v 2 2 

= 0) = d (^(-3 + 6z - 3z 2 )(K 2 ) 2 - (1 - 2zf (1 3) 2 [1 3] 2 + (6 - 18* + \2z 2 )K 2 [1 3] (1 3) 
For q = 1: 

(1 - 2z) 2 [3 1] (1 3) - (1 - 2z)K 2 



—B mi) (s) 



s=0 2 (3 4) 2 [12] 2 (l-2z) 

= C x ( (1 - 2zf [1 3] 2 (1 3) 2 + (-1 + 7 -z - 3z 2 )K 2 (1 3) [1 3] 



(6.70) 



d »(2;1;1) 



to 



Ci ( (4 - Wz + 6z 2 )(i^ 2 ) 2 + 2(1 - 2zf [1 3] 2 (1 3) 2 + (-10 + 30z - 2\z 2 )K 2 [1 3] (1 3)' 



s=0 
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For q = 2: 
1 d 2 



2ds 2 * 2 >° {S) 



rz((l-2z) [3 1]<1 3) -(l-z)K 2 ) 



(6.71) 



=o (3 4) 2 [1 2] 2 

= Ci f(-l + 2z - z 2 )(K 2 ) 2 - (1 - 2z) 2 (1 3) 2 [1 3] 2 + (4 - Uz + \2z 2 )K 2 (1 3) [1 3] 



s=0 



All together, we get the following result for the N\ term. 

5(A' 2 ) 4 



clK]m _ mmw: (i 

[121 2 (4 3) 2 V2 



+ 



[1 2r(4 3) 

K 2 [l 3] (1 3) / 1 2 



4z(l - ^ 



(4 3) 2 [1 2] 2 V 6 3 



z (l — z ) J ^ i 

3[1 2] 2 (4 3) 2 6[1 2] 2 (4 3) 2 

( - - , ( 1 : ) ) a 3) 2 [1 3] 2 + (~ - 3*(1 - z))K 2 (1 3) [1 3] 



(6.72) 



The result of cr 12 |34] 
Final bubble coefficient: 

2 f(l 3)[1 3] (if 2 ) 3 (I 

c [12|34] ^ 



(^ 2 ) 2 \ [12] 2 (4 3) 2 \2 

K 2 [l 3] (1 3) f.l 2 



4z(l - z) 



3 [1 2] 2 (4 3) 2 



6[1 2] 2 (4 3) 2 



+- 



Q - ~-z{\ - z)) (1 3) 2 [1 3] 2 + Q- 3z(l - z))K 2 (1 3) [1 3] 



1 



+ 



(4 3) 2 [1 2] 2 V'6 3 
(12) [3 4] [4 3] 



K 2 [1 2] 



(1 2) 2 [3 4] 2 



-4(m 2 + n 2 )s + 6(m 2 + fi 2 )t - 2st) 



(6.73) 



3s 2 t 2 

where we used the definitions K 2 = t, (1 3) [1 3] = — s — t. 

6.5 Comparison with the literature 

The i-channel cut of A(l~ , 2~~, 3 + , 4 + ) admits a decomposition in terms of cuts of master integrals as 
shown in Fig|| Its expression was given in Eq.(5.33) of [26|, and reads 



A| ermion (l-,2-,3 + ,4- 



t— cut 



-2 A| caJaJ (l-,2-,3 + ,4 H 



t— cut 



(4vr 



.2-e 



AT e (tJ A -I 2 (t)) 



(6.74) 



t— cut 



with 



and 



^| calar (l-,2-,3 + ,4 + ) 



t-cut (47T 



1 ^(lip>*~~ + itu) 



t— cut 



(1 2)< 



.(1 2) 2 [3 4] 2 



(1 2) (2 3) (3 4) (4 1) 



K 2 s 



12 



(6.75) 



(6.76) 
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where we neglected the cut-free term, I\ and .^(O). In standard notation we have s = S12, t = S23, u - 
-s-t = S13. 

Now we translate the expression of [^] into our canonical basis, using the identities of Appendix |A|. 

1 



4 ermion (l-,2-,3 + ,4 H 



2^ calar (l~,2~,3 + ,4 + ) 



t— cut 



^ ree (^4-/ 2 (t)) 



t— cut 



(4 



t-cut (47r) 2 - £ 



.77) 



(6.78) 

t— cut 



,; .IV" ( - ( |/ 2 BFM [1] + -h ~ -/ 2 BFM [m 2 + ft 2 ]) + -I 2 BFM [m 2 + M 2 ] 



--/ 4 BFM [(m 2 + ^i 2 ) 2 ] + t/ BFM [m 2 + /j 2 ] - / 2 BFM [1] 



(6.79) 



t— cut 



(1 2) 2 [3 4] 2 (2 



st 



6 St s 



+ ^/BFM [(m 2 + M 2 )2] _ t/ BFM [m 2 + ^ 



(6.80) 



t— cut 



We have reproduced every one of these coefficients, up to an overall minus sign in the amplitude. 

7. From polynomials in u to final coefficients 

As proven in the Appendix ^, the coefficients of 2-, 3-, and 4-point functions in four dimensions 
are polynomials in u (or equivalently fi 2 ), of known degree d: for boxes, d = \{n + 2)/2]; for triangles, 
d = [(n + l)/2]; for bubbles, d = [n/2]; where [x] denotes the greatest integer less than or equal to x. Using 
this fact, we can generally represent any coefficient of the master integral as 



(7.1) 



r=0 



The coefficients c r are in one-to-one correspondence to the coefficients of the shifted- dimension master 
integrals (see Section 2). 

To compute the c r analytically, one can proceed with the standard differentiations with respect to u, at 
u = 0: 

1 d r 



k\ du' 1 



-Pd(u) 



(7.2) 



u=Q 



When the differentiations are time consuming, or the analytic expression is not needed, one can switch to 
the following numerical procedure, and extract the c r algebraically, by projections. 

1. Generate the values Pd,k, (k = 0, d — 1), 



Pd,k = Pd(uk) , 



(7.3) 
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by evaluating Pd(u) at particular points: 



(7.4) 



2. Using the orthogonality relations for plane waves, one can obtain the coefficient c r simply by the 
following formula: 



cr = ~£ P d , k e 2mrk/d . (7.5) 



k=0 
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A. Change of Basis 



To compare our results to the literature, we need to convert the master integrals used in [26, 27] to our 
canonical (4 — 2e)-dimensional basis, ( |2.1| ). For clarity, we now denote the basis used in this paper by 
I BFM , while the other integrals in this appendix are defined according to |26|, p7fl . The first point is then 
that 



I n = i(-l)"(47rr- e C M (A.l) 
We use the identities from Appendix A. 4 of [26] to perform the conversion. 

J 4 = / 4 [ m 2 + ^] = i(47r) 2 - e / 4 BFM [m 2 + fj 2 ] (A.2) 
I 2 (t) = i(4vr) 2 - e / 2 BFM [l] (A.3) 

I^ DM = '-hit) + \h ~ ~ 3 4 h3) = .(4vr) 2 - (^/ 2 BFM [1] + \h - \ir U [m 2 + /, 2 ]) (A.4) 

j(l,3) = J (l,3) [m 2 + At 2 ] = . (47r) 2- ei BFM [m 2 + M2] (A5) 

K 4 = h[(m 2 + /i 2 ) 2 ] = i(4vr) 2 -/ 4 BFM [(m 2 + /i 2 ) 2 ] (A.6) 
B. The ^-dependence of the coefficients 

Here we analyze the u-dependence of the integral coefficients given by our formulas. First, we prove that 
they are polynomials in u. Then, we present some alternate formulas where this polynomial dependence 



is more explicit. This material is a straightforward generalization of the analysis in the massless case [30 
so we omit many of the details here. 
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To begin, we rewrite our vectors Rj, Qj from ( 0.23 ) and ( 3.24| ) in the following way: 

r2 



Rj = ~ ((1 " 2z) + Pi + ft*, 



P^[P3-^K). (B.I) 



+ (B . 2) 

/ M 2 -M 2 \ 

Qi = -((!- 2z) + ^ 2 2 j <& + a,-*, (B.3) 

(JT,- - if) 7 ^ 2 *' 2x ^ 2 1 ™ 2 



K 2 



/ M 2 -M 2 \ Kj + Mf-mj 



Notice that 



Pj-K = 0, qj-K = 0. (B.6) 



Using (gig ), we have 



E,-(u) = -p(y/T=u) Pj + f3jK, Qj(u) = -P(VT^) qj + aj K. (B.7) 

This is the same expression as in the massless case, except for the factor (5. The point is that now all 
it-dependence is in the factor yl — u, just as in the massless case. In fact, we should now consider the 
factor f3y/l — u as our basic quantity. The proof that the integral coefficients are polynomials in u was 
performed by considering the (demonstrably finite) series expansion in y/1 — u, and showing that the odd 
powers drop out. Therefore, the same arguments now carry over to the series expansion in f3y/l — u. 

B.l Triangle coefficients 

Let us begin with triangle coefficients. The null vectors P S) j exhibit a simple dependence on u. Specifically, 

PsM = -(3(VT^H)P qa;i , 

where 



P qs4 = q s ±\ \I^t\K, (B.8) 

which is manifestly independent of u. In defining the spinor components of P s j, we can place the "in- 
dependent factor inside the antiholomorphic spinor, i.e., 



\P s ,i) = \P q J , \P Si \ = -P(VT=u)\P q ,j]. (B.9) 
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Then, for the triangle coefficients, we have 
{K 2 ) 1+n 1 



C[Q S ,K] 



2 

d n+l 

rU + l 



-/3^r^)-+i(v /Z 4e^)™+ 1 (n + 1)! (P qsA P qs , 2 ) n+1 
Ujti (P gs ,i - rP qat2 \R,(u)Q s (u)\P gsil - rP qs>2 



dT n+i \j^ =x ^ s {PqsA _ T P qs , 2 \Q t ( U )Q s ( U )\P qsA - rP qst2 
Further, we make use of some identities, 

(i\Qt(u)Q s (u)\£) = U{Qt(u) - ^Q s (u))Q s (u)\£) = -py/T^Ufa - ^q s )Q s ( u )\£ 
(£\R t (u)Q s (u)\£) = -PVT^(£\(jh - £q.)Q a (u)\£\ 
along with the definitions 

qt = (qt q s ), Pt = [Pt q s )- 

a s a s 

Our final form for the triangle coefficient is: 

n+l 



(B.10) 



T=0 



(B.ll) 



C[Q S ,K] = ~ 



K 2 
=4^ 



(n + l)!(P.i P s , 2 y 



n+l 



d n+l I U^liPsA-rPsMPjQsHlPs,!-^ 



8,2 



dT n+l \j[tl-#s ( P sA ~ TP.,2\qtQ.(v.)\P.,l ~ tP., 2 ) 



+ {P a>1 <-> P s>2 } 



(B.12) 



r=0 



Here, the u-dependence is concentrated entirely within the vector Q s (u), since we have made sure to choose 
the spinor components wisely in ( |B.9| ), so that the holomorphic spinors are u- independent. 

B.2 Bubble coefficients 



We follow the same procedure as with triangles, and make use of the same definitions ( [B.8D , (B.9). The 
■u-dependence can be concentrated within the vectors Rj(u) and Q r (u), Q p (u), along with the explicit factor 
of \Jl — u, in the following formulas. 



1 d« 



k n 



C[K] = {K 2 )^ £ B^ n _ q (s) + £ £ {Kn=r(s) - B^lf\s] 

q=0 \ r=l a=q 



(B.13) 



s=0 



where 



d n 



i (2 V .K)^ YES m(*){K - *n)\4 



dr« I n\[ri\nK\n\« (t + (l ^"+1 ft* (£\Q p ( u )(K - srj)\£) 



\£)^\K-ri)\ v l 



,(B.14) 



r=0 
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B 



(r;6;l) 



n .t 



(*) = 



1) 



6+1 



t+1 



1 (P r ,i - TPrMPg^ 

t+1 



b l { _P^T—^ )b+ i v /_ Aq 2 K 2 b+1 (P p )^rf ^(t+ 1) (P r>1 - rP ri2 |K|P gril ] 



(P r ,i - rP r , 2 |Q r (u)77|P r ,i - rP r , 2 ) b n"=i fc (fi.i ~ rP r , 2 |^(n)(^ - sr/j^i - rP^) 
(P r ,l - rP r , 2 |T7Jr|P r ,i - rP r / +1 rip=i,^r <^,i - TP r , 2 \Q p (u)(K - srf)\P r>l - rP r , 2 ) 



,(B.15) 



r=0 



<f ;2) (^) 



1) 



6+1 



t+1 



1 <J ? r,2-rP r ,i|7/|P gr , 2 ] 

t+1 



6! (_ / VT^)6+i v /_ 4g 2^2 b + 1 ( p rji p r2) 6 dr b ^(i + 1) (P r . 2 - rP ril |i^|P gri2 



(P r<2 - rP r ,i|Q r ( M )7y|P r , 2 - rP r ,i) fe n"=i fc (fi,2 - rP^i^)^ - S 7y)|P r , 2 - rP r ,i) 
X (Pr,2 - rP r!l \ V K\P ri2 - rP r ,i) n+1 Up=i, P ^r ( p r,2 ~ tP t ,i\Q p (u)(K - s V )\P r , 2 - rP r ,i) 

B.3 Box and pentagon coefficients 



.(B.16) 



T=0 



Although the formula ( [4.1| ) for box and pentagon coefficients looks simple, the u-dependence now gets 
complicated. We consider the separate cases k = 2, k = 3, and k > 4. 

B.3.1 The case k = 2 

In this case, there is only one box, and no pentagons. The box coefficient is given by 



(K 



2\2+n 



rij=l \ P{Qj (u),Qj(«));l ( u ) I Rj ( u ) I P(Qj {u),Qi(u));2 («) 

71+2 + l P (O i (tt) I O i («)) i l(«) "(Qj(u),Qi(u));2( ? W-4 ' J 



V \- p (Qi(«),Q<(«));i( w )l- K 'l- p (Oi(«),<3i(«))j2(«) 



Given the vectors Qj, Qj,K that select a particular box, it is useful to construct a vector 
is orthogonal to all three, and independent of u: 



(qi,qj,K) 



that 



if 2 



(B.18) 
(B.19) 



As in the massless case, the u-dependence can be concentrated in a single factor, aS gi ' q i\u). If all input 
quantities are set to their values with u = 0, except for adjusting the definition of R s (u) as follows, 



R s (u) 



R,(u) 



(gi,9j,^)\ 2 



/3 2 (l-u) + 



4if 2 [a i a J (2 gi - gj )-aj ! g J-a|g a 2 ] 
(2< ?l -<7,) 2 -4g l 2 g, 2 



2 i IK^Xonajl^qvq^-cx^-o^qW 



(B.20) 



(B.21) 



/? 2 + 



(2 9i . % )2-4gf g ! 



then the value of the box coefficient remains the same. 
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In summary, the box coefficient for k = 2 is given by 

rn+2 

i L'- i-- " / I 



{K 2 )2+ n mii{P j iil\R.(u)\P jii2 



(P ji;1 \K\P ji;: 



in+2 



+ {P ji;1 «-> P,;. 2 } 



(B.22) 



where 



and 



In evaluating ( B.23 ), it is useful to observe the following: 



( 

VPs ■ q 



(q t ,q 3 ,K), _ e(p s ,qi,qj,K) _ e(P s , Kj, Kj, K) 



K 2 



K 2 



(B.23) 



(B.24) 



(B.25) 



This formula ( B.22 ) looks the same as in the massless case; the difference is the appearance of (3 in ( |B.23| ), 
both explicitly and through the definition ( B.21| ) of ec-H'ii) . 



B.3.2 The case k = 3 

Here there is a pentagon, as well as three boxes. The differences from the massless case are all based in 
the definitions of Rj (u) , Qj (u) : there is always a factor of accompanying v'l — u, and mass parameters 
enter into the definitions (BJ2), (B.5) of (3j,Uj. 



When we make these adjustments, we find that the pentagon coefficient takes the same form as in the 
massless case, 



n+3 



but the definition of (3. 



C[Q l ,Q 3 ,Q t ] = (K 2 f +n \[^ qM) 

s=l 

now includes mass parameters: 



(B.26) 



D{qi>qj>qt;ps) _ a{Ki,Kj,Kt;P a ) 



(B.27) 



pvm,qj,qt,Psj _ 

(K 2 + M 2 - m 2 )e(P s ,K„K, K t ) + (K 2 + M 2 - m 2 )e(K h P S ,K, K t ) 

K 2 e(Ki, Kj,K, K t ) 
(g2 + M 2 _ M%)e(Kj, Kj , P s , K t ) + (K 2 + M 2 - m 2 )e(Kj, Kj, K, P 8 ) 

K 2 e(Ki, Kj,K, K t ) 

The expression ( |B.27 ) is symmetric in Ki, Kj, Kt, K (recall that Mi is the mass associated with K in this 
context). 
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The box coefficients are given by 

{K 2 ) 2+n 

C[Qi, Qj]k=3 = 5 

n+3 

-<n« 

s=l 



(q,, qj ,q t ; Pa ) (Pji;l\K \Pjip) 



P ii;1 |Q t (n)|P ii;2 



+ {P ji;1 <-» P ii;2 } 



(B.28) 



The derivation of ( |B-2§| ) involved the result from the case k = 2. All mass-dependence is already included 
in the definitions ( |B,21| ) and ( |B.23| ), along with the similarly defined vector Qt(u): 

R s (u) 



Qt(u) 



^^(a(«*)(u) - l)(-^r^) + R s (u = 0) 



B.3.3 The case k > 4 

In the derivation of the formulas, we introduce the following functions: 

uK^KM (Kf + Ml - m 2 )e(K, Kj,K s ,K t ) + (K 2 + M 2 - mj)e(K u K, K s ,K t 



(B.29) 
(B.30) 



7s 



K 2 e{Ki,K h K,K t ) 
(K 2 + Mf - m2)e(K i , Kj , K, K t ) + (K 2 + M 2 - m 2 )e(K u Kj,K s , K) 



+ 



K 2 e(Ki, Kj,K, K t ) 



e{K u K 3 ,K,K t ) 



(B.31) 



The numerator of ^ s Ki,K ^ Ka ' Kt ^> \ s symmetric in Ki, Kj, K s , Kt; the denominator breaks this symmetry by 
singling out K s . 

We find the following results. 

The pentagon coefficients are given by 



C^Q^Qt] = (K ) 

The box coefficients are given by 

{K 2 ) 2+n 



T-rn+3 n(qi,qj,qr,Ps) 
2\3+n Lls=l Ps 



i-rfc (Ki,Kj,K w ,Kt) ' 

1 lw=l,w^i,j,t iu> 



(B.32) 



n"i n (^ii^(«)i^ ;2 



C[Qi, Qj]k>4 



E 



(P i i ; l|^|P i d n+2 nt Wi j(^ ; l|4MI^ ; 2 



ji;2\ 



* — ' TT^ (Ki,Kj;K m ,K t ) I D \r\ I Wry 

-\Pji\X *-* Pji;2S 



(B.33) 
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Again, all the u-dependence is concentrated in R(u) and Q(u). The definitions of Pji- a , Rs(u), Qt(u), 
P^'if'P°\ and ^K f ,K w ,K t ) g . ven . n ft£2$ t (|R29|) , (gjg ), ( ET27| ), and (EOlD , respectively. 
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